Genetic determinants of SARS-CoV-2 and the clinical outcome of COVID-19 in Southern Bangladesh.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic has had a severe impact on population health. The genetic determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in southern Bangladesh are not well understood. METHODS: This study aimed to determine the genomic variation in SARS-CoV-2 genomes that have evolved over 2 years of the pandemic in southern Bangladesh and their association with disease outcomes and virulence of this virus. We investigated demographic variables, disease outcomes of COVID-19 patients and genomic features of SARS-CoV-2. RESULTS: We observed that the disease severity was significantly higher in adults (85.3%) than in children (14.7%), because the expression of angiotensin-converting enzyme-2 (ACE-2) diminishes with ageing that causes differences in innate and adaptive immunity. The clade GK (n = 66) was remarkable between June 2021 and January 2022. Because of the mutation burden, another clade, GRA started a newly separated clustering in December 2021. The burden was significantly higher in GRA (1.5-fold) highlighted in mild symptoms of COVID-19 patients than in other clades (GH, GK, and GR). Mutations were accumulated mainly in S (22.15 mutations per segment) and ORF1ab segments. Missense (67.5%) and synonymous (18.31%) mutations were highly noticed in adult patients with mild cases rather than severe cases, especially in ORF1ab segments. Moreover, we observed many unique mutations in S protein in mild cases compared to severe, and homology modeling revealed that those might cause more folding in the protein's alpha helix and beta sheets. CONCLUSION: Our study identifies some risk factors such as age comorbidities (diabetes, hypertension, and renal disease) that are associated with severe COVID-19, providing valuable insight regarding prioritizing vaccination for high-risk individuals and allocating health care and resources. The findings of this work outlined the knowledge and mutational basis of SARS-CoV-2 for the next treatment steps. Further studies are needed to confirm the effects of structural and functional proteins of SARS-CoV-2 in detail for monitoring the emergence of new variants in future.

The structural basis of hyperpromiscuity in a core combinatorial network of type II toxin-antitoxin and related phage defense systems.

Toxin-antitoxin (TA) systems are a large group of small genetic modules found in prokaryotes and their mobile genetic elements. Type II TAs are encoded as bicistronic (two-gene) operons that encode two proteins: a toxin and a neutralizing antitoxin. Using our tool NetFlax (standing for Network-FlaGs for toxins and antitoxins), we have performed a large-scale bioinformatic analysis of proteinaceous TAs, revealing interconnected clusters constituting a core network of TA-like gene pairs. To understand the structural basis of toxin neutralization by antitoxins, we have predicted the structures of 3,419 complexes with AlphaFold2. Together with mutagenesis and functional assays, our structural predictions provide insights into the neutralizing mechanism of the hyperpromiscuous Panacea antitoxin domain. In antitoxins composed of standalone Panacea, the domain mediates direct toxin neutralization, while in multidomain antitoxins the neutralization is mediated by other domains, such as PAD1, Phd-C, and ZFD. We hypothesize that Panacea acts as a sensor that regulates TA activation. We have experimentally validated 16 NetFlax TA systems and used domain annotations and metabolic labeling assays to predict their potential mechanisms of toxicity (such as membrane disruption, and inhibition of cell division or protein synthesis) as well as biological functions (such as antiphage defense). We have validated the antiphage activity of a RosmerTA system encoded by Gordonia phage Kita, and used fluorescence microscopy to confirm its predicted membrane-depolarizing activity. The interactive version of the NetFlax TA network that includes structural predictions can be accessed at http://netflax.webflags.se/.

Context-based sensing of orthosomycin antibiotics by the translating ribosome.

Orthosomycin antibiotics inhibit protein synthesis by binding to the large ribosomal subunit in the tRNA accommodation corridor, which is traversed by incoming aminoacyl-tRNAs. Structural and biochemical studies suggested that orthosomycins block accommodation of any aminoacyl-tRNAs in the ribosomal A-site. However, the mode of action of orthosomycins in vivo remained unknown. Here, by carrying out genome-wide analysis of antibiotic action in bacterial cells, we discovered that orthosomycins primarily inhibit the ribosomes engaged in translation of specific amino acid sequences. Our results reveal that the predominant sites of orthosomycin-induced translation arrest are defined by the nature of the incoming aminoacyl-tRNA and likely by the identity of the two C-terminal amino acid residues of the nascent protein. We show that nature exploits this antibiotic-sensing mechanism for directing programmed ribosome stalling within the regulatory open reading frame, which may control expression of an orthosomycin-resistance gene in a variety of bacterial species.

A hyperpromiscuous antitoxin protein domain for the neutralization of diverse toxin domains.

Toxin-antitoxin (TA) gene pairs are ubiquitous in microbial chromosomal genomes and plasmids as well as temperate bacteriophages. They act as regulatory switches, with the toxin limiting the growth of bacteria and archaea by compromising diverse essential cellular targets and the antitoxin counteracting the toxic effect. To uncover previously uncharted TA diversity across microbes and bacteriophages, we analyzed the conservation of genomic neighborhoods using our computational tool FlaGs (for flanking genes), which allows high-throughput detection of TA-like operons. Focusing on the widespread but poorly experimentally characterized antitoxin domain DUF4065, our in silico analyses indicated that DUF4065-containing proteins serve as broadly distributed antitoxin components in putative TA-like operons with dozens of different toxic domains with multiple different folds. Given the versatility of DUF4065, we have named the domain Panacea (and proteins containing the domain, PanA) after the Greek goddess of universal remedy. We have experimentally validated nine PanA-neutralized TA pairs. While the majority of validated PanA-neutralized toxins act as translation inhibitors or membrane disruptors, a putative nucleotide cyclase toxin from a Burkholderia prophage compromises transcription and translation as well as inducing RelA-dependent accumulation of the nucleotide alarmone (p)ppGpp. We find that Panacea-containing antitoxins form a complex with their diverse cognate toxins, characteristic of the direct neutralization mechanisms employed by Type II TA systems. Finally, through directed evolution, we have selected PanA variants that can neutralize noncognate TA toxins, thus experimentally demonstrating the evolutionary plasticity of this hyperpromiscuous antitoxin domain.

Sal-type ABC-F proteins: intrinsic and common mediators of pleuromutilin resistance by target protection in staphylococci.

The first member of the pleuromutilin (PLM) class suitable for systemic antibacterial chemotherapy in humans recently entered clinical use, underscoring the need to better understand mechanisms of PLM resistance in disease-causing bacterial genera. Of the proteins reported to mediate PLM resistance in staphylococci, the least-well studied to date is Sal(A), a putative ABC-F NTPase that-by analogy to other proteins of this type-may act to protect the ribosome from PLMs. Here, we establish the importance of Sal proteins as a common source of PLM resistance across multiple species of staphylococci. Sal(A) is revealed as but one member of a larger group of Sal-type ABC-F proteins that vary considerably in their ability to mediate resistance to PLMs and other antibiotics. We find that specific sal genes are intrinsic to particular staphylococcal species, and show that this gene family is likely ancestral to the genus Staphylococcus. Finally, we solve the cryo-EM structure of a representative Sal-type protein (Sal(B)) in complex with the staphylococcal 70S ribosome, revealing that Sal-type proteins bind into the E site to mediate target protection, likely by displacing PLMs and other antibiotics via an allosteric mechanism.

FlaGs and webFlaGs: discovering novel biology through the analysis of gene neighbourhood conservation.

SUMMARY: Analysis of conservation of gene neighbourhoods over different evolutionary levels is important for understanding operon and gene cluster evolution, and predicting functional associations. Our tool FlaGs (standing for Flanking Genes) takes a list of NCBI protein accessions as input, clusters neighbourhood-encoded proteins into homologous groups using sensitive sequence searching, and outputs a graphical visualization of the gene neighbourhood and its conservation, along with a phylogenetic tree annotated with flanking gene conservation. FlaGs has demonstrated utility for molecular evolutionary analysis, having uncovered a new toxin-antitoxin system in prokaryotes and bacteriophages. The web tool version of FlaGs (webFlaGs) can optionally include a BLASTP search against a reduced RefSeq database to generate an input accession list and analyse neighbourhood conservation within the same run. AVAILABILITY AND IMPLEMENTATION: FlaGs can be downloaded from https://github.com/GCA-VH-lab/FlaGs or run online at http://www.webflags.se/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A holin/peptidoglycan hydrolase-dependent protein secretion system.

Gram-negative bacteria have evolved numerous pathways to secrete proteins across their complex cell envelopes. Here, we describe a protein secretion system that uses a holin membrane protein in tandem with a cell wall-editing enzyme to mediate the secretion of substrate proteins from the periplasm to the cell exterior. The identity of the cell wall-editing enzymes involved was found to vary across biological systems. For instance, the chitinase secretion pathway of Serratia marcescens uses an endopeptidase to facilitate secretion, whereas the secretion of Typhoid toxin in Salmonella enterica serovar Typhi relies on a muramidase. Various families of holins are also predicted to be involved. Genomic analysis indicates that this pathway is conserved and implicated in the secretion of hydrolytic enzymes and toxins for a range of bacteria. The pairing of holins from different families with various types of peptidoglycan hydrolases suggests that this secretion pathway evolved multiple times. We suggest that the complementary bodies of evidence presented is sufficient to propose that the pathway be named the Type 10 Secretion System (TXSS).

A widespread toxin-antitoxin system exploiting growth control via alarmone signaling.

Under stressful conditions, bacterial RelA-SpoT Homolog (RSH) enzymes synthesize the alarmone (p)ppGpp, a nucleotide second messenger. (p)ppGpp rewires bacterial transcription and metabolism to cope with stress, and, at high concentrations, inhibits the process of protein synthesis and bacterial growth to save and redirect resources until conditions improve. Single-domain small alarmone synthetases (SASs) are RSH family members that contain the (p)ppGpp synthesis (SYNTH) domain, but lack the hydrolysis (HD) domain and regulatory C-terminal domains of the long RSHs such as Rel, RelA, and SpoT. We asked whether analysis of the genomic context of SASs can indicate possible functional roles. Indeed, multiple SAS subfamilies are encoded in widespread conserved bicistronic operon architectures that are reminiscent of those typically seen in toxin-antitoxin (TA) operons. We have validated five of these SASs as being toxic (toxSASs), with neutralization by the protein products of six neighboring antitoxin genes. The toxicity of Cellulomonas marina toxSAS FaRel is mediated by the accumulation of alarmones ppGpp and ppApp, and an associated depletion of cellular guanosine triphosphate and adenosine triphosphate pools, and is counteracted by its HD domain-containing antitoxin. Thus, the ToxSAS-antiToxSAS system with its multiple different antitoxins exemplifies how ancient nucleotide-based signaling mechanisms can be repurposed as TA modules during evolution, potentially multiple times independently.

ABCF ATPases Involved in Protein Synthesis, Ribosome Assembly and Antibiotic Resistance: Structural and Functional Diversification across the Tree of Life.

Within the larger ABC superfamily of ATPases, ABCF family members eEF3 in Saccharomyces cerevisiae and EttA in Escherichia coli have been found to function as ribosomal translation factors. Several other ABCFs including biochemically characterized VgaA, LsaA and MsrE confer resistance to antibiotics that target the peptidyl transferase center and exit tunnel of the ribosome. However, the diversity of ABCF subfamilies, the relationships among subfamilies and the evolution of antibiotic resistance (ARE) factors from other ABCFs have not been explored. To address this, we analyzed the presence of ABCFs and their domain architectures in 4505 genomes across the tree of life. We find 45 distinct subfamilies of ABCFs that are widespread across bacterial and eukaryotic phyla, suggesting that they were present in the last common ancestor of both. Surprisingly, currently known ARE ABCFs are not confined to a distinct lineage of the ABCF family tree, suggesting that ARE can readily evolve from other ABCF functions. Our data suggest that there are a number of previously unidentified ARE ABCFs in antibiotic producers and important human pathogens. We also find that ATPase-deficient mutants of all four E. coli ABCFs (EttA, YbiT, YheS and Uup) inhibit protein synthesis, indicative of their ribosomal function, and demonstrate a genetic interaction of ABCFs Uup and YheS with translational GTPase BipA involved in assembly of the 50S ribosome subunit. Finally, we show that the ribosome-binding resistance factor VmlR from Bacillus subtilis is localized to the cytoplasm, ruling out a role in antibiotic efflux.

In silico identification and characterization of common epitope-based peptide vaccine for Nipah and Hendra viruses.

OBJECTIVE: To explore a common B- and T-cell epitope-based vaccine that can elicit an immune response against encephalitis causing genus Henipaviruses, Hendra virus (HeV) and Nipah virus (NiV). METHODS: Membrane proteins F, G and M of HeV and NiV were retrieved from the protein database and subjected to different bioinformatics tools to predict antigenic B-cell epitopes. Best B-cell epitopes were then analyzed to predict their T-cell antigenic potentiality. Antigenic B- and T-cell epitopes that shared maximum identity with HeV and NiV were selected. Stability of the selected epitopes was predicted. Finally, the selected epitopes were subjected to molecular docking simulation with HLA-DR to confirm their antigenic potentiality in silico. RESULTS: One epitope from G proteins, one from M proteins and none from F proteins were selected based on their antigenic potentiality. The epitope from the G proteins was stable whereas that from M was unstable. The M-epitope was made stable by adding flanking dipeptides. The 15-mer G-epitope (VDPLRVQWRNNSVIS) showed at least 66% identity with all NiV and HeV G protein sequences, while the 15-mer M-epitope (GKLEFRRNNAIAFKG) with the dipeptide flanking residues showed 73% identity with all NiV and HeV M protein sequences available in the database. Molecular docking simulation with most frequent MHC class-II (MHC II) and class-I (MHC I) molecules showed that these epitopes could bind within HLA binding grooves to elicit an immune response. CONCLUSIONS: Data in our present study revealed the notion that the epitopes from G and M proteins might be the target for peptide-based subunit vaccine design against HeV and NiV. However, the biochemical analysis is necessary to experimentally validate the interaction of epitopes individually with the MHC molecules through elucidation of immunity induction.